Formation of molecular oxygen in ultracold O + OH reaction 
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We discuss the formation of molecular oxygen in ultracold collisions between hydroxyl radicals 
and atomic oxygen. A time-independent quantum formalism based on hyperspherical coordinates is 
employed for the calculations. Elastic, inelastic and reactive cross sections as well as the vibrational 
and rotational populations of the product O2 molecules are reported. A J-shifting approximation 
is used to compute the rate coefficients. At temperatures T = 10 — 100 mK for which the OH 
molecules have been cooled and trapped experimentally, the elastic and reactive rate coefficients 
are of comparable magnitude, while at colder temperatures, T < 1 mK, the formation of molecular 
oxygen becomes the dominant pathway. The validity of a classical capture model to describe cold 
collisions of OH and O is also discussed. While very good agreement is found between classical 
and quantum results at T = 0.3 K, at higher temperatures, the quantum calculations predict a 
larger rate coefficient than the classical model, in agreement with experimental data for the O + 
OH reaction. The zero-temperature limiting value of the rate coefffcient is predicted to be about 



6 X 10" 



cm molecule 



a value comparable to that of barrierless alkali-metal atom - dimer 



systems and about a factor of five larger than that of the tunneling dominated F -I- H2 reaction. 
PACS numbers: 



I. INTRODUCTION 



Important experimental progress is being made in cre- 
ating ultracold molecules in tightly bound vibrational 
levels Very recently, formation of ground state 

molecules in the vibrational level v ~ has been re- 
ported by different groups for homonuclcar molecules 
such as Cs2 01 , Rb2 Q , as well as for heteronuclear polar 
molecules such as RbCs [i], KRb and LiCs 0. There 
has also been much progress in the measurement of rate 
coefficients of barrierless reactions involving alkali-metal 
atoms at cold and ultracold temperatures. This includes 
atom - molecule collisions such as Rb -I- Rb2 [a, B, Cs 
+ Cs2 [13, [Hi, Na + Na2 [11, Rb/Cs + RbCsllI and 
molecule - molecule collisions such as Cs? + Cs9 |lOl.fl^. 
Na2 + Na2 [11, and Rba -t- Rba [J. The typical 
rate coefficients of these reactions are on the order of 
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-11 



10- 



s depending on the col- 



lisional system, the temperature and the vibrational lev- 
els probed. All these experimental studies indicate that 
inelastic and reactive processes occur at significant rates 
at ultracold temperatures, in accordance with theoretical 
predictions on barrier reactions [3, [H, [13, [H, [H, H^l 
as well as a number of alkali-metal atom - dimer sys- 
tems [m, m m [n [n m, [13, [n, [Si such as li -i- 

Li2, Na -I- Na2 and K -I- K2 which proceeds without an 
energy barrier in the entrance channel. The alkali-metal 
systems arc characterized by triatomic complexes with 
deep potential wells which make them challenging sys- 
tems for accurate quantum calculations. For these heavy 
systems, the density of triatomic states is large and it 
leads to strong couplings between them, enhancing in- 
elastic quenching or reactive scattering [sij . However, 



explicit quantum calculations of molecule - molecule sys- 
tems involving alkali-metal atoms are computationally 
intractable. In a recent theoretical study of vibrational 
relaxation in collisions between II2 molecules, Qucmcncr 
et al. [3I showed that the relaxation rate coefficients at 
ultralow temperature can attain large values for some 
near-resonant processes which involve simultaneous con- 
servation of internal energy and total internal rotational 
angular momentum of the colliding molecules. 

While there have been a number of theoretical stud- 
ies of ultracold reactive collisions of tunneling dominated 
reactions with chemically distinct reactants and prod- 
ucts [H, [H, [13, [H, [H, [2I, there have been no such 
studies involving barrierless chemical reactions at ultra- 
cold temperatures. Here we investigate the exothermic 
reaction 



0(3p) + 0H(2h) ^ H(2S) + Oaf 



(1) 



at cold and ultracold temperatures as an example of a 
barrierless chemical reaction involving non-alkali-metal 
atom systems. The reaction is of key interest in oxygen 
chemistry in the interstellar medium, OH chemistry in 
the upper stratosphere and mesosphere, and combustion 
chemistry (see Ref. [H and references therein). The OH 
radical has also been cooled and trapped using the buffer 
gas and stark decelerator techniques [H, [H, [H, [stI . 
High precision measurements of its radiative lifetime [38[ 
and its hyperfine constant [s^, [11 have recently been 
reported. An experimental study of scattering between 
cold OH molecules and He atoms and D2 molecules has 
recently been reported [4l|. The cooling and trapping 
studies have stimulated a number of theoretical investi- 
gations involving the OH molecule in the last few years. 
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Gonzalez-Sanchez et al. [4^ reported rotational relax- 
ation and spin-flipping of OH in He collisions at ul- 
tralow energies. They found that rotational relaxation 
occurs more efficiently than elastic collisions at vanish- 
ing collision energies. By carrying out quantum calcu- 
lations of Rb -I- OH collisions at ultracold temperatures 
Lara et al. [i^, [31 explored the possibility of sympa- 
thetic cooling of OH molecules by collisions with Rb 
atoms. They argued that efficient sympathetic cool- 
ing of OH molecules in the ground vibrational state by 
collisions with Rb atoms is unlikely to occur due to 
the large inelastic rate coefficient. External fields can 
also have important effects on ultracold molecular colli- 
sions [i^, [H, [ij . Avdeenkov and Bohn investigated the 
effect of external fields on ultracold collisions between OH 
or OD molecules ^,^1131. Ticknor and Bohn [5l| also 
studied OH— OH collisions in the presence of a magnetic 
field. They showed that magnetic fields of several thou- 
sand Gauss suppress inelastic collision rates by about 
two orders of magnitude. In a recent study [s^, we re- 
ported quantum dynamics calculations of reaction ([1]) for 
T = 10 — 600 K and found no significant decrease of the 
rate coefficient in the temperature range 39— 10 K, in ac- 
cordance with conclusions of a recent experimental study 
by Carty et al. [s^l. Our calculations for reaction prob- 
abilities were in excellent agreement with those of Xu et 
al. [5^ for collision energies Ec > 0.012 eV. 

In this paper, we present the quantum dynamics of 
reaction ([1]) at low and ultralow collision energies to ex- 
plore the behavior of complex forming chemical reactions 
at cold and ultracold temperatures. Since OH molecules 
have been experimentally cooled and trapped at low tem- 
peratures, we believe that collisional properties of the O 
-I- OH reaction will be of considerable interest to the cool- 
ing and trapping community. The paper is organized as 
follows: The details of the quantum mechanical formal- 
ism along with convergence tests are discussed in Section 
II. Results of cross sections, rate coefficients, and state- 
to-state product distributions are given in Section III. We 
also include in this section a discussion on the usefulness 
of a classical model in describing cold collisions of O and 
OH. Conclusions are presented in Section IV. 



II. METHODOLOGY 

A. Potential energy surfaces 

We employed a modified version of the electronically 
adiabatic ground state (1 '^A") potential energy surface 
(PES) of HO2 calculated by Kendrick and Pack ^ 
using a diatomics-in-moleculc (DIM) formalism. This 
new version includes improvements to the long-range be- 
havior and is referred to as the DIMKP PES. In par- 
ticular, a switching function was implemented which 
more smoothly "turns-on" the long-range van der Waals 
potentials for the electronic ground states of both O2 
and OH. This switching function is given by /switch = 



0.5(tanh(Q!(r — rp)) -l- 1) where a = 1, tq = 7.0 aq for 
02(^Sg ) and tq = 10 oq for 0H(2H). A minor global re- 
fitting of the DIM HO2 PES to the original set of ab initio 
data was required in order to account for the new switch- 
ing functions and ensure a smooth transition to the long- 
range behavior. The same long-range coefficients were 
used as in the original version of the surface, for OH(^H): 
CgO-H = 9.295 Ehal Ct"" = 169.09 E^al CO-" - 
4060.85 Ehal", and for 02(^2"): C^"° = 14.89 Eha^ 
CgO-o = 206.67 Ehal, and CPq-^ = 3753.745 £:,,aj° [H]. 
The global fit in the interaction region is essentially iden- 
tical to the original fit with nearly the same rms devia- 
tion of 0.099 eV (2.3 kcal/mol). The improvements to the 
long-range behavior are important for the ultracold col- 
lisions studied in this work but should not significantly 
affect the results of previous scattering calculations at 
higher (thermal) energies [11]. We present in Fig. [T] the 
potential energy curves for the 3 lowest ^A" states of 
HO2 for the DIMKP PES, for a 0-HO Uncar approach. 
We also note in the inset of Fig. [1] that the DIMKP PES 
predicts a shallow conical intersection along the 0-HO 
approach due to the crossing of the OH(H) and OH(E) 
states. For a fixed ron = 1-83 Qq, this crossing occurs at 
i?o_Ho = 10.8235 ao and its energy Hes 8.30 x 10""* eV 
(« 9.6 K) below the asymptotic energy of the O + OH 
channel. For comparison purposes, we also employed the 
ab initio PES computed by Xu, Xie, Zhang, Lin, and 
Guo [13, [Hi, referred to as the XXZLG PES. The XX- 
ZLG PES has been used in a number of quantum dynam- 
ics calculations of the O -I- OH system at high collision 
energies [H, [H, [H . The present DIMKP PES was em- 
ployed for the first time for this reaction in our previous 
work [ssj and it is preferred at low energies as it includes 
accurate long-range coefficients. 
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FIG. 1: Potential energy curves for the 3 lowest ^A" states of 
HO2 for the DIMKP PES, for a O-HO linear approach with 
roH — 1.83 ao. The conical intersections arising at 7?o-ho ~ 
10.8 ao is shown in the inset. 
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B. Quantum mechanical approach and convergence 

tests 

The quantum dynamics calculations have been per- 
formed using the adiabatically adjusting principle-axis 
hyperspherical (APH) approach of Pack and Parker [b^] • 
The method uses the democratic Smith- Whittcn hyper- 
spherical coordinates in the inner region (p < 17 ao) that 
includes the triatomic well of the HO2 system and the 
Delves-Fock hyperspherical coordinates in the outer re- 
gion (p > 17 ao) in the valleys of the H -I- O2 and O -I- 
OH arrangement channels. For a given value of the to- 
tal angular momentum quantum number, J, and hyper- 
spherical radius, p, the wavefunction is expanded onto 
a basis set of adiabatic functions, which are eigenfunc- 
tions of a triatomic hyperangular Hamiltonian. A hy- 
brid DVR/FBR primitive basis set [6l[ combined with an 
Implicity Restarted Lanczos algorithm is used to diago- 
nalize the hyperangular Hamiltonian matrix. The time- 
independent Schrodinger equation yields a set of differen- 
tial close-coupling equations in p, which are solved using 
the log-derivative matrix propagation method of John- 
son [62|- The log-derivative matrix is propagated to a 
matching distance (pm) where asymptotic boundary con- 
ditions are applied to evaluate the reactance matrix 
and the scattering matrix S'-' . The square elements of the 
S'' matrix provide the state-to-state transition probabili- 
ties, P'' . The matching distance and all other parameters 
employed in the calculations were determined by opti- 
mization and extensive convergence studies. To secure 
convergence of the reaction probabilities, the number of 
hyperspherical channels, n, included in the close coupling 
equations is 393. This is sufficient to obtain converged 
results at low and ultralow energies. The cross sections 
are calculated using the standard formulae: 



Pm (ao) 26.1 



32.7 



39.9 



44.0 



-10 
-9 



(eV) 
10" 
10- 
10"* 
10-^ 
10-*^ 
10"^ 
10"* 
10"^ 



0.4995 
0.4965 
0.4877 
0.4596 
0.3708 
0.1432 

0.005536 0.005214 0.005300 0.005217 
0.001531 0.001516 0.001506 0.001507 



0.5019 


0.4786 


0.4679 


0.4977 


0.4742 


0.4635 


0.4853 


0.4610 


0.4505 


0.4472 


0.4206 


0.4110 


0.3369 


0.3085 


0.3021 


0.1164 


0.1108 


0.1123 



TABLE I: Convergence of the elastic cross section cr/i " in 
units of 10-i^cm^ with the matching distance p,„ for O + 
OH(i; = 0, j = 0) using the DIMKP PES, for n = 393. 



Pm (ao) 


26.8 


32.7 


39.9 


44.0 


(eV) 


cr/e"° X v 








10-10 


0.1965 


0.2777 


0.3123 


0.3133 


10-9 


0.1954 


0.2755 


0.3095 


0.3104 


10"* 


0.1920 


0.2688 


0.3010 


0.3018 


10-^ 


0.1819 


0.2490 


0.2760 


0.2765 


10"'' 


0.1544 


0.1979 


0.2123 


0.2119 


lO"'^ 


0.09835 


0.1080 


0.1080 


0.1073 


10-* 


0.03748 


0.03713 


0.03724 


0.03714 


10"=* 


0.008704 


0.008745 


0.008761 


0.008765 



TABLE II: Convergence of the non-thermal reactive rate co- 
efficients, a:^^'^ X V, in units of 10^^" cm"^ molecule-^ s'^ with 
the matching distance pm for O -I- OH(ii : 
DIMKP PES, for n = 393. 



0, j = 0) using the 



= -^11 -s 



J\2 



k^ 4" 

reactive j 



J |2 



where i, j denote initial and final quantum states. The 
non-thermal elastic and reactive rate coefficients are 
given by CToi x v and a-^c x v where v = ?ik/p is the 
incident velocity for relative motion of the O atom and 
the OH molecule. 

The convergence of the elastic cross section cr^i^" and 
the non-thermal reactive rate coefficient a'^^^ x v are 
presented in Tab. U and Tab. |TT] for O + OYi{v = 0, j = 
0) on the DIMKP PES. At vanishing collision energies, 
these quantities attain finite values as required by the 
Bethe-Wigner laws [glliil: 



const. 



X V 



const. 



The tables also show convergence of the results with the 
matching distances. At low energies the elastic cross sec- 
tion converges less rapidly with the matching distance 
Pm than the reactive one. This is because the long-range 
contribution to the interaction potential is not negligi- 
ble compared to the kinetic energy in the entrance chan- 
nel, even for moderately large values of the hypcrradius. 



Thus, elastic cross sections are generally more sensitive 
to the long-range tail of the interaction potential. For 
the parameters in Tab. |T] and Tab. [H] the reactive rate 
coefficients are converged to within 1% while the elastic 
cross sections are converged within 3%. All of the final 
results presented here are obtained using p,„ = 44.0 ao, 
which is especially necessary to get converged results at 
< 10-4 eV Hi. 



III. RESULTS AND DISCUSSION 

A. Cross sections and rate coefficients 

The J = elastic and reactive cross sections are plot- 
ted as a function of the collision energy in Fig. [2] for the 
DIMKP PES. The corresponding non-thermal rate co- 
efficients are presented in Fig. [31 In both figures, the 
Bethe-Wigner laws are satisfied at ultralow collision en- 
ergies. The elastic cross section and the reactive rate co- 
efficient converge to values of 4.6 x lO"^"' cm^ molecule"^ 
and 3.1 x 10"^^ cm"^ molecule--'^ s"^ respectively in the 
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— O + OH(v=0,i=0) -> O + OH(v=0,i=0) 
O + OH(v=0,i=0) ^ O + OH(v=0,Z j*0) 

— O + OH(v=0,i=0) ^ H + Oj 
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FIG. 2: (Color online) J = cross section as a function of 
the collision energy for the elastic (red curve), inelastic (green 
dashed curve) and reactive (black curve) collisions of O atoms 
with ORiv = 0, j = 0) molecules for the DIMKP PES. 




FIG. 3: (Color online) J = non-thermal rate coefficient as 
a function of the collision energy for the elastic (red curve), 
inelastic (green dashed curve) and reactive (black curve) col- 
lisions of O atoms with OH(t; = 0,j = 0) molecules for the 
DIMKP PES. 



FIG. 4: (Color online) Same as Fig. |3] but for the XXZLG 
PES. 



O + OH(v=0,j=0) ^ H + 0^{v,) ; £^=10"'' K 
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12 3 

FIG. 5: "Vibrational populations of the product O2 molecule 
in O 4- OH(ii = 0,j = 0) H + 02(1^/) reaction evaluated 
using the DIMKP PES at a collision energy of 10"*^ K and 
J = 0. 



limit of zero energy. The clastic rate coefficient and the 
reactive cross section behave respectively as the square 
root and the inverse of the square root of the colli- 
sion energy. The Bethe-Wigner regime is reached at 
Ec « 10~^ K, where the reactive rate coefficient is of 
about one order of magnitude higher than that of the 
elastic counterpart. At Ec = 10~® K the reactive rate 
coefficient is two orders of magnitude larger than the 
elastic one. This is reminiscent of other exothermic bar- 
rierless systems such as alkali-metal atom - diatom col- 
lisions [III, m, m, [13, M, M, At energies 
between Ec = 0.01 K and 30 K, elastic and reactive rate 
coefficients are of comparable magnitude. For Ec > 
the elastic scattering is more efficient than reactive scat- 
tering. The inelastic rotational excitation to the j = 1 
level opens up at i?c = 52 K and its cross section is com- 
parable to the reactive contribution. 



The J = contribution to the elastic, inelastic and re- 
active rate coefficients evaluated using the XXZLG PES 
are presented in Fig.[4]as functions of the collision energy. 
For i?c > 1 K, all three rate coefficients are of compara- 
ble magnitudes with those obtained from the DIMKP 
PES. For Ec < 1 K, the reactive rate coefficient is about 
two orders of magnitude smaller than that obtained us- 
ing the DIMKP PES. The differences can be traced to 
the incomplete description of the long-range interaction 
potential in the XXZLG PES. It docs not properly in- 
clude the long-range potential in the O + OH channel as 
it was not designed for quantum dynamics calculations 
at ultralow energies. Thus, the comparison of the results 
obtained using the two surfaces highlights the crucial role 
of the long-range interaction potential in chemical reac- 
tion dynamics at low temperatures. The importance of 
the long-range intermolecular forces in the 0-OH system 
in determining its rate coefficient has been pointed out 
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O + OH(v=0,j=0) ^ H + 0„(v,,j,) ; E =10 K 
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FIG. 6; Rotational populations of the product O2 molecule 
in O + OH(u = 0,j = 0) H + 02(u/ = 0, 1,2,3, j/) reac- 
tion evaluated using the DIMKP PES at a collision energy of 
10'*' K and J = 0. 



by Clary and Werner [66|. In the following, we restrict 
our calculations to the DIMKP PES. 

The vibrational populations of molecular oxygen 
formed in O + OH(u = 0,j ~ 0) collisions, evaluated 
using the DIMKP PES at an energy of 10"*^ K and 
J = are presented in Fig. [5l Oxygen molecules are 
predominantly formed in vibrational levels Vf = [0—2] 
with a slight preference for the ground vibrational state 
Vf = 0. Thus, a significant fraction of the O2 molecules 
formed is in excited vibrational states. The rotational 
level distributions of molecular oxygen formed in O + 
OH(v = 0,j = 0) collisions for J = are presented in 
Fig- El for each final vibrational level Vf populated by the 
reaction at an energy of 10~^ K. The results show that 
low rotational levels are generally preferred for each of 
the final vibrational levels. Since the incident channel 
includes only the s-wave, the final rotational distribution 
is largely determined by the anisotropy of the interaction 
potential. 

Explicit calculations of rate coefficients would require 
reaction probabilities for all contributing values of the 
total angular momentum quantum number (J). This is 
a computationally demanding problem for the O + OH 
reaction if full quantum dynamics calculations are em- 



ployed, especially when many J values contribute to the 
reaction probability. The J-shifting approximation [67j is 
widely used to compute rate coefficients when full quan- 
tum calculations are not practical. This is a good ap- 
proximation for barrier reactions which involve a transi- 
tion state but not for complex forming reactions. Nev- 
ertheless, the J-shifting approximation has been applied 
to the O -|- OH reaction in a number of previous stud- 
ies [H, [H, m, [gl, [7^ and it has been demonstrated 
that it can predict rate coefficients within about 40% 
of numerically exact calculations [s^]. Here we use the 
J-shifting approximation [gJ to compute the rate coef- 
ficients for the O -I- OH reaction. The rate coefficient is 
given by the expression: 



X ( V(2J + l)e-<-/('=-^)^ 
X / P::f="{Ec) e'^^/^"-^^ dE, (2) 



where ks is the Boltzmann constant, Pl'^f ° is the J ^ 







reaction probability and E^YiUt i^ the height of the effec- 
tive barrier for a given partial wave J in the entrance 
channel. To determine the barrier height for a given par- 
tial wave, we first evaluate the effective potential, T/^: 



off 



;iV(j + 1) 

2/i(i?0-OH)2 



+ Knin(-Ro-OH) 



(3) 



where VJnin(fio-OH) is the minimum energy path of the 
reaction as a function of the atom - molecule center-of- 
mass separation, i?o-OH, and ^ is the reduced mass of 
the 0-OH system. In Eq. ([2]), Qr = Qtrans x Qoi is the 
reactant partition function. For the translational par- 
tition function we used the standard formula, Qtrans = 

^^^2 J . For the electronic partition function we used 
the expression given by Graff and Wagner [7l|: 



(5 + 3e-228/T ^ e-326/r)(2 + 2e-205/^) 



cl — 



The effective barriers, V^^, for the DIMKP PES are 
shown in Fig. [7| for J = [0 - 200; 10]. As Fig. [7] il- 
lustrates, the "reef visible in the effective potential at 
Ro-on = 5 ao for low values of J becomes an effective 
barrier for J > 70 as indicated by the bold line. The bar- 
rier heights £'s'^ift and their locations i?fhift are reported 
in table Im] for J = [1 - 10]. 

The elastic, inelastic, and reactive rate coefficients of 
the O + OH(w = 0, j = 0) reaction evaluated using 
the J-shifting method for the DIMKP PES are shown 
in Fig. El for T = 10"'' „ iq3 k. As indicated earlier, the 
Bethe-Wigner regime is reached for temperatures below 
T « i?c « 10~^ K. The value of the rate coefficient in 
the zero-energy limit is 6.2 x 10~^^ cm'^ molecule"^ s~^. 
This is a factor of 5 smaller than the rate coefficient 



reported as cr, 



j=o 



V in Fig. [31 The difference comes 
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T(K) 



FIG. 7: Effective potentials V^g as functions of -Ro-oh for J = 
[0 - 200; 10] for the DIMKP PES. The bold line corresponds 
to J = 70. 



FIG. 8: (Color online) Rate coefficients of O + On{v = 0, j = 
0) for the elastic (red curve), inelastic (green dashed curve) 
and reactive (black curve) collisions for the DIMKP PES. The 
rate coefficient from a classical Langevin model (black dashed 
curve) is also shown (see text for detail). 



from the electronic partition function in the denomi- 
nator in Eq. ^ which is equal to 5 as T — > 0. In 
experiments using the stark decelerator methods, OH 
molecules were cooled to T = 10 — 100 mK. In this 
temperature range, our computed values of the elas- 
tic rate coefBcients are comparable to the reactive ones 
and will not favor sympathetic cooling of OH by col- 
lisions with O atoms. Similar conclusions have been 
found by Lara et al. [i^, [3] for collisions between OH 
molecules with Rb atoms. The relatively large rate co- 
efficient for the reaction in the zero-temperature limit 
indicates that barrierless exothermic reactions occur at 
significant rates at ultracold temperatures, in agreement 
with similar results for alkali-metal atom - diatom reac- 
tions [m, [13, m, [13, [H, [H, [13, [H, [H. Fig. E shows 

that the minimum value of the rate coefficient for the O 
-|- OH reaction is about 4.9 x 10~^^ cm"^ molecule"^ s^^ 
atT = 5xlO~^K. This provides a lower limit for O2 
formation by reaction llj . We note that the inelastic pro- 
cess becomes more probable than the reactive process for 
T > 330 K. 



B. Classical capture model 

We shall now compare our quantum dynamics re- 
sults with a classical capture model also known as the 
Langevin model [t^] . This model has been shown to work 
quite well for a certain range of collision energy, for atom- 
exchange reactions such as K -t- K2 [1^, Li -I- Li2 24 1, 
and also for non- reactive Rb + OH collisions [H, 44 1. 



We have shown [3l[ that the Langevin model predicts 
rate coefficients in close agreement with experimental val- 
ues for Rb + Rh2 i, Cs -t- Cs2 [13, [HI, and Rb/Cs + 
RbCs [1^ collisions reported recently. For the O + OH 
reaction. Clary and Werner [6^ used an adiabatic cap- 
ture theory [73| , while Davidsson and Nyman 7^ used a 
generalized Langevin model. Both studies showed good 



agreement between theoretical and experimental rate co- 
efficients. Here we test the validity of the Langevin model 
against a TIQM and J-shifting method for the O + OH 
reaction. For a non-rotating diatomic molecule (j = 0) 
at large atom - molecule separations -Rq-oh, for a given 
value of J, the interaction potential can be approximated 
by an effective potential, ^La^g, composed of a repulsive 
centrifugal term and an attractive van dcr Waals term: 



* Lang 



ftV(J+l) CeO-OH 

2^(i?0-OH)2 (i?0-OH)' 



At an atom - diatom distance of 



-"Lang 



OH 



r9 J{J+1) 



1/4 



the height of the effective potential is given by: 



-'^Lang 



3^/3 



(C°-°H)- 



1/2 



J(j + i) 

2fl 



3/2 



The quantities E^^^^ and R'l^^^g are reported in Tab. IIIII 
for J = I - 10. The (7^"°" coefficient is evaluated by 
fitting the long-range part of the minimum energy path 
V^ni„(i?o-OH). This yielded a value of C^-°" « 9.2 Eha^ 
which is much smaller than the sum of the atom - atom 
coefficients C^-^ + C^-° = 24.2 Eha^. In the DIM 
model, the smaller C^"*-*^ coefficient is due to the pres- 
ence of significant diatomic mixing in the OH diatomic 
states for roH = L83 ao (see for example Eqs. 34 — 37 in 
Ref. [IHl). This mixing gives rise to an effective Ce coef- 
ficient of the 0-HO approach of ai C^-^ + as C^"° = 
9.2 EhttQ where the multiplicative constants ai and 02 
depend on the level of diatomic mixing. As roH increases, 
the diatomic mixing decreases and ai — > 1 and a2 ^■ 
The dependence of the long-range coefficients for atom - 
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diatom interactions on the diatomic separation has also 
been noted and investigated in recent work on the Li 
+ Li2 and Na + Na2 systems [t^ [t^. For the atom - 
molecule reduced mass we used fi = 15023.74 au. In the 
Langevin model, the rate coefficient as a function of the 
temperature is given by the expression: 



J 


ELft (K) R 


shift (^o) 


ELns (K) R 


Lang (S-o) 


1 


0.022 


25.4 


0.022 


25.4 


2 


0.11 


19.3 


0.11 


19.3 


3 


0.30 


18.6 


0.32 


16.3 


4 


0.54 


18.4 


0.68 


14.3 


5 


0.86 


18.0 


1.26 


12.9 


6 


1.26 


17.5 


2.08 


11.9 


7 


1.75 


17.1 


3.20 


11.1 


8 


2.34 


16.7 


4.68 


10.4 


9 


3.04 


16.4 


6.53 


9.8 


10 


3.84 


16.1 


8.83 


9.3 



TABLE III: Comparison of the heights and positions of the 
effective barriers between the J-shifting approximation and 
the Langevin model. 

The rate coefficient predicted by the Langevin model 
is plotted in Fig. [5] (black dashed curve). Except in the 
Bethe-Wigner regime for T < 10~^ K, where the classi- 
cal model is not valid, the Langevin model predicts rate 
coefficients in semi-quantitative agreement with those 
obtained from the quantum calculations. Overall, the 
reactive rate coefficients oscillate slightly around the 
Langevin line, as in other barrierless systems such as K 
+ K2 dH and Li -I- Li2 H^]. However, the Langevin 
model yields good agreement only for a restricted range 
of temperatures. The lower limit of the model prediction 
is restricted by the number of partial waves included in 
the calculations. In previous studies [H, [23i it has been 
found that when three or more partial waves are included, 
the quantum and classical results are in good agreement. 
If less than three partial waves are involved, the quantum 
character becomes dominant and the results cannot be 
compared with the classical model. These previous stud- 
ies have also shown that the maximum of the quenching 
rate coefficient for a given partial wave J occurs at about 
a collision energy comparable to the height of the barrier, 
^Lang- Thus, the lower limit corresponds to a collision 
energy of El=^^ = 0.32 K (see TabUm for the present 
system. The upper limit is bounded by the long-range 
part of the potential. When the Langevin radius i^L^ng 
located at a distance where the long-range part of the po- 
tential is not described by the van der Waals interaction, 
^iang ''^ill differ from -E^jft . In this case, the classical re- 
sults will differ from the quantum calculations. Tab. IIIII 



shows that this is the case for J = 3 and the upper limit 
of the model also corresponds to a collision energy of 
^iang = 0-32 K. Thus, the Lan gevin model gives quan- 
titative agreement for temperatures around T « 0.3 K, 
as seen in Fig. [S] For T > 0.3 K, the rate coefficient 
depends on the exact details of the effective potential in 
the entrance channel. However, the classical result is in 
overall agreement with the quantum result, though the 
classical model predicts a smaller value for the rate co- 
efficient. This is because E^^^^^^ < E^^^^ for J > 3 (see 
Tab, imp and reactivity is less probable to occur with the 
Langevin model. Therefore, the Langevin model provides 
a lower limit of the rate coefficients for T = 1 — 100 K 
which affirms the theoretical conclusions of Ref . [s^ and 
experimental conclusions of Ref. [s^ , that the rate coef- 
ficient of reaction ^ is unlikely to vanish for T < 10 K. 
We also note that the Langevin rate coefficient is in semi- 
quantitative agreement with the inelastic rate coefficient 
in the range T = 100—1000 K. The difference is less than 
30 %. 



IV. CONCLUSION 

This paper presents the first quantum mechanical in- 
vestigation of an ultracold barrierless reaction with chem- 
ically distinct reactants and products. We investigated 
dynamics of molecular oxygen formation in ultracold col- 
lisions of the hydroxyl radical and atomic oxygen using 
a time-independent quantum mechanical method based 
on hyperspherical coordinates. It has been found that 
formation of molecular oxygen occurs with a relatively 
large rate coefficient of 6.2 x 10~^^ cm^ molecule"^ s~^ 
at ultracold temperatures. The oxygen molecules are 
mainly formed in the v = — 2 vibrational levels with 
a slight preference for the v = level. Calculations 
show that at temperatures of T = 10 — 100 mK, the 
elastic cross sections are not large enough to achieve 
efficient evaporative cooling in collisions between OH 
molecules and O atoms. We predict a lower limit of 
4.9 X 10~^^ cm'^ molecule"^ s~^ for the rate coefficient for 
the O + OH(i; = 0, j = 0) reaction at T = 5 x lO'^ K. 
This shows that formation of O2 molecules is significant 
even at ultracold temperatures. It has been found that a 
classical capture model is valid for temperatures around 
T « 0.3 K where quantum and classical calculations yield 
comparable results. Based on our analysis of the long- 
range interaction potential we find that for T = 1 — 10 K 
the Langevin model can provide a lower limit to the 
quantum reactive rate coefficients calculated within the 
J-shifting approximation. 

Future work will consider the effects of the geomet- 
ric phase and non-adiabatic couplings between different 
PESs of HO2 on the reaction dynamics. The geometric 
phase due to the conical intersection of two PESs may 
have an important effect at low and ultralow collision 
energies where only a few partial waves contribute to the 
reaction probabilities. 
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